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It is shown that globally-coupled oscillators with pulse 
interaction can synchronize under broader conditions than 
widely believed from a theorem of MiroUo & Strogatz Q. 
This behavior is stable against frozen disorder. Beside the 
relevance to biology, it is argued that synchronization in relax- 
ation oscillator models is related to Self-Organized Criticality 
in Stick-Slip- like models. (LPTHE preprint: 94/29) 

PACS numbers: 05.45+b,Q5.40.-|-j 

Large assemblies of oscillator units can spontaneously 
evolve to a state of large scale organization. Collective 
synchronization is the best known phenomenon of this 
kind, where after some transient regime, a coherent os- 
cillatory activity of the set of oscillators emerges. This 
effect has attracted much interest in biology for the study 
of large scale rhythms in populations of interacting ele- 
ments The south-eastern fireflies, where thousands of 
individuals gathered on trees fiash altogether in unison, 
is the most cited example |^]. Other examples are the 
rhythmic activity of cells of the heart pacemaker, of cells 
of the pancreas, and of neural networks . Most of the 
works on synchronization have used models in which the 
interaction between the oscillators is smooth and contin- 
uous in time (see e.g. |^). Comparatively, few analytical 
results are known on models where the interactions are 
episodic and pulse-like although they are rele- 

vant to several biological situations as fireflies and neural 
networks. This letter deals with the emergence of syn- 
chronization in a very simple model of globally coupled 
Integrate and Fire (IF) oscillators, see Eq.(^ for the def- 
inition and It is shown that synchronization occurs 
under broad conditions on the properties of the oscilla- 
tors, thus generalizing a theorem of MiroUo and Strogatz 
[Q, the usual interpretation of which restricted the syn- 
chronization conditions in a too drastic way. Further- 
more I show that the synchronized state is stable against 
a frozen disorder of the oscillator properties, a subject 
much studied in models with continuous interaction pi. 
This result is different from that of a previous study W] 
on a closely related model. 

Beside synchronization another form of collective or- 
ganized behavior is known to occur in large assemblies 
of elements with pulse interaction, that is Self-Organized 
Criticality (SOC). This concept has been proposed in ||] 
to describe out of equilibrium systems that are gener- 



ically critical, i.e. that organize into a scale invariant 
critical state spontaneously, i.e. without tuning of a con- 
trol parameter. Systems displaying SOC are externally 
driven with a drive slower than any other characteristic 
time. These models are made critical by the choice of a 
threshold dynamics that forbids them to follow adiabat- 
ically the external drive. They can be modelled as cou- 
pled map lattices, and can be divided into two subclasses 
based on the concept of oscillators. In the subclass (a) the 
external drive acts globally and continuously on all the 
lattice sites of the coupled map, until one of them reaches 
the threshold, in which case it relaxes to zero: each site 
is therefore a relaxation oscillator. These are the stick- 
slip- like models [|| . They are deterministically |^ or not 
|lO| , p| driven systems, for which no conserved quantity 
is known, and can exhibit SOC even when they are dis- 
sipative The systems of the subclass (b) are driven 
at each time step by the increment of a unique site, so 
that sites are not oscillators. These are the sandpile-like 
models HJl^ . The model studied in this letter is a simple 
model of fireflies flashing and is a mean field version of 
class (a) models. As this mean-field approach suggests, 
we argue in the following that there is a close relationship 
between SOC and synchronization pO| , pl] , p^ . 

The biological mechanisms controlling the fireflies 
flashing have been extensively studied (see j2j for a re- 
view) by testing the response of an isolated insect to sin- 
gle or periodic flashes of light. It appears that there are 
several mechanisms for synchronization and entrainment 
among the different species, so that no general conclusion 
can be drawn about both the mechanism and the bio- 
logical function of synchronized flashing. It is therefore 
interesting in this context to study under which general 
conditions a synchronization may emerge in relaxation 
oscillator models. Experiments indicate that the rhyth- 
mic spontaneous flashing of the firefiies, is in general gov- 
erned by a flash-control pacemaker in the brain, cycling 
as an oscillator with a given intrinsic rhythm between 
a basal level of excitation (state variable of the oscilla- 
tor) and a fully excited flash-triggering level After 
each flash the state variable resets quickly to the basal 
level. Light-flash stimuli on a single firefiy have differ- 
ent effects depending on the species. This has motivated 
the introduction of essentially two models, the phase ad- 
vance and phase delay models In the phase advance 
model, the effect of a pulse of light is to advance the 
flashing of each insect pushing the state variable towards 
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the threshold. Flashing delay relative to the intrinsic 
rhythm is impossible in this case. This accounts for the 
behavior of at least two species: Photinus pyralis and 
Photinus concisus^ which however display only transient 
synchrony possibly because these species are rover [^. 
On other species, as P. cribellata, a pulse of light inhibits 
the oscillator, thus delaying the next flash. This is inter- 
preted in the phase delay model as the resetting of the 
state variable to the basal level. In this model, advance 
or delay of flashing is possible. Let us mention that other 
fireflies species have behaviors that cannot be described 
by any of these two models. The model we study in the 
following ||l|,|l^ is equivalent to the phase advance model. 
The mean field model. It consists of N relaxation oscil- 
lators Oi represented when they evolve freely by a state 
variable Ei — E{t) € [0, Ec = 1] monotonically increas- 
ing in time with period 1. The interaction is so defined 
that when Ei > 1 it relaxes to zero and increments all 
the other oscillators by a pulse a/N: 

where a £ [0, 1] is a dissipation parameter. The model 
is taken into the slow drive limit by assuming that any 
avalanche of firings, i.e. any succession of firings triggered 
by the relaxation of one oscillator, is instantaneous. We 
assume a supplementary rule that in most circumstances 
does not change the behavior of the system but simpli- 
fies the discussion. This is the so-called absorption rule 
of 0, which is equivalent to a reaction death time of 
the oscillators immediately after their relaxation. This is 
taken into account by assuming that oscillators whithin 
the same avalanche are not incremented by the pulses 
resulting from the following relaxations in the avalanche. 

The theorem of MiroUo & Strogatz 0, that applies 
to this model, states that for a convex function Ei — 
E{t) the system synchronizes completely but for a set 
of initial conditions of Lebesgue measure null. Although 
the validity of this theorem cannot be questioned, we will 
see in the following that for all the physically interesting 
situations, the system synchronizes even for a linear and 
a range of concave functions E{t). 

Let us call configuration the set of ordered distinct val- 
ues e[^^ < E^^ < . .. < Em] = 1 of the state variables 
present in the system just before the (j -I- l)-th avalanche. 
To each E^ corresponds a group of oscillators at 

this value (X]i='i ^i"'^ = Let call cycle, the time nec- 
essary for all the rrij groups to avalanche exactly once. 
To trace the evolution of the system, it is useful to fol- 
low, cycle after cycle, the gaps sf ^ — £'|:i_\— E^p'' between 

a) 

the values of successive groups. If one of these gaps ' 
becomes smaller than the value aN^'i^-^/N of the pulse 
of the (i -f- l)-th group, then the (i)-th group gets ab- 
sorbed by the {i + l)-th group. For a linear E{t), an ele- 
mentary calculation shows that the gap between a group 



[i + 1) and a smaller one (i) is reduced during one cycle 
by 5 Si = \a{Ni — Nij^i)/N\. Therefore large groups un- 
avoidably absorb the smaller ones that follow them and 
become larger and larger it is this positive feedback 
that leads to synchronization. The only way for stopping 
this mechanism and therefore the evolution towards syn- 
chronization is to get a configuration where all groups 
are of equal size. Apart in the random initial conditions, 
where all the groups are of size one, it is exceptional for 
large N to be stuck during the evolution in such a state. 
It is in fact a difficult task to calculate in general the 
probability for this situation to occur, but, first, exten- 
sive numerical calculations show that this is indeed ex- 
ceptional even for small N , second, it is even a physically 
ill-defined question since it depends on the factorization 
of iV, and third, if one imagines that the oscillators are 
not exactly identical, the slightest splitting between their 
pulse intensities forbids any 5si to vanish. Therefore to 
see if the system does not get stuck in the initial config- 
uration, the only relevant question is to know if it is pos- 
sible with random initial conditions to get groups larger 
than one during the first cycle. The probability for pair 
formation may easily be computed from the Poissonian 
distribution P{s)ds = Ne~^'ds, {N > 1), of gaps s be- 
tween the N random numbers E'l^^ in [0, 1]. The number 
of oscillator pairs in the random initial values satisfying 
the synchronizing conditions s < a/N is then: 

N / P{s)ds = N{1 - e""), AT > 1 (2) 
Jo 

The absorption process of oscillators into groups takes 
place as long as there is at least one such a gap. It fol- 
lows directly from (Q) that this is typically the case when 
a > 1/N [|o|. Initial configurations of values where no 
gap is smaller than a/N are in principle possible, but for 
large systems with a 0(e^"^) probability of occurence. 
In conclusion, the set of linear oscillators completely syn- 
chronizes almost always, when a is of order 1 compared 
to N. The reason why this is not in contradiction with 
the theorem of MiroUo & Strogatz, is that their recursive 
demonstration requires that two single oscillators syn- 
chronize. While this is the case for convex oscillators, 
it is not so for linear and concave ones. However, I dis- 
agree with the ususal interpretation of the theorem of 
MiroUo & Strogatz as the impossibility of synchroniza- 
tion for non convex oscillators: in fact, this theorem tells 
nothing in this case. In Q, it was noted that convex- 
ity is not necessary for synchronization in the case of a 
model with transmission delays. We see here that this 
condition is not even necessary. Moreover, even concave 
oscillators can synchronize, provided that the concavity 
is not too large. In this case, positive feedback is in com- 
petition with concavity, that has for effect to increase 
the gap between groups as they approach the threshold. 
For small concavities positive feedback prevails, whereas 
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for larger concavity, groups are not able to grow. This 
process can be seen on the class of concave functions 
E{t) ^t'',a> 1. Let us call = 1 - (1 - na/iV)(i/'^) 
the value of the phase difference between two successive 
groups of n and m oscillators (n > m) below which the 
first group absorbs the second one. The first return map 
of the phase difference Xi between these two groups when 
the largest one gets to the threshold (i.e. the state value 
of the second group is (1 — Xi)"") has an attractive fixed 
point xq. The condition for synchronization is xo < Xc- 
For a fixed pulse a/N there is an upper value for the con- 
cavity amax = 1 — log(l — a/iV)/ log 2 below which the 
preceding condition is fulfilled for all group sizes n > m. 
In this case the system synchronizes as with linear oscilla- 
tors. Since Umax IS the most stringent bound, we expect 
that synchronization almost always occurs for a larger 
range of values of a. This has been extensively verified 
numerically [pT| . 

Up to now, all the oscillators were identical. However 
it is known that firefly populations have a spread of indi- 
vidual frequencies Therefore it is necessary to show 
that the phase delay model accounts for synchronization 
even in this case. 

Frequency distribution. Different frequencies may be al- 
lowed in IF models, by keeping an identical threshold 
for the oscillators, but letting them have different slopes 
Ei{t) — Git , Ei G [0, 1], so that the a^'s are the internal 
frequencies. Positive feedback is here also the effective 
mechanism leading to synchronization for small disorder 
of frequencies. For larger frequency dispersion the for- 
mation of large groups is not always possible, since two 
oscillators relaxing together at one moment do not nec- 
essarily relax again together at the next cycle. In fact 
two questions must be answered: first, can stable groups 
exist, and second, can such groups be formed during the 
evolution? As we will see, this second point is not prob- 
lematic, and stability is the most stringent condition for 
the emergence of synchronization. Let us first notice 
that a necessary condition for the stability of a group 
of n oscillators, is that the relaxation of the fastest one 
triggers the avalanche of all the others. Therefore sta- 
ble groups must fire at the rhythm of the elements with 
highest frequency. Let us now study the stability of a 
group of (n+1) elements in the background of the other 
— (n -t- 1) oscillators, assuming that the subset of the 
n first ones, Oi, . . . , On of frequencies ai > . . . > a„ is 
stable, and that the (n + l)-th oscillator 0„+i, has the 
lowest frequency Un+i- Assuming that Oi, . . . , 0„, On+i 
have just relaxed together, the condition for again relax- 
ing together after a cycle is: 
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where A is the summed strength of the pulses of all the 
other oscillators. Therefore the most stringent condition 
of stability for a group of n -|- 1 oscillators is 



^^^<(i-l)-^ Vz = l...n+I (4) 

oi A 

Therefore, given a system of A^ oscillators with a random 
uniform distribution of frequencies centered on a with 
width D — Umax ~ o,mim the probability that the whole 
system remains stable is: 
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where p = N/D and S — aamax/N. One can see on 
Fig. 1^ that even for large D/a there is a relatively high 
probability to have a realization of frequencies allowing 
complete synchronization as long as a is not too small, in 
which case the system decouples and synchronization is 
less probable. One can verify that, when the probability 
Pn is not vanishingly small, it is almost independent of 
A^ for large A^ (~ 100). Now we will see that any two 
stable groups of oscillators will, at some moment, relax 
together. Although positive feedback is still acting as 
in the case of identical oscillators, it is in competition 
with disorder that tends to desynchronize groups. How- 
ever, contrary to the case of identical linear oscillators, 
it is easy to show that the dispersion of frequencies im- 
plies that the relative positions of any two groups change 
monotically in time. Therefore at some moment the two 
groups must relax together. If their union is itself stable, 
it forms a unique stable group. This proves that if the 
distribution of frequencies fulfills the stability conditions, 
Eq.(^, the system completely synchronizes. If this con- 
dition is not satisfied, there is partial synchronization: 
the subsystems that individually fulfill Eq.(|^) synchro- 
nize. Extensive numerical simulations confirm all these 
conclusions 

It is natural to think of other kinds of frozen disorder. 
For instance, the pulses do not necessarily need to be 
of the same strength, and also the oscillators may have 
different threshold values. As expected, and for similar 
reasons as those previously explained for the case of a 
spread of frequencies, synchronization occurs also in this 
case and even for strong disorder. These results are dif- 
ferent from the conclusions in [Q where for a similar but 
different model, the completely synchronized state is un- 
stable even against weak disorder. 

Lattice models. Beside synchronization or quasi- 
synchronization, lattice models of pulse-coupled oscilla- 
tors may display SOC. This is the case for the models of 
class (a), for which however SOC seems to be related to 
a tendency to synchronization |^,^. SOC also appears 
when a system is perturbed, which otherwise should syn- 
chronize totally or partially |10 1^, or which should be 



3 



periodic [|5j. The most striking example of class (a) 
models is the Olami-Feder-Christensen (OFC) model ^ 
that consists of oscillators Ei on a square lattice, that re- 
lax to zero when they exceed a given threshold i^c, thus 
incrementing their nearest neighbors by a pulse which is 
a {a < 1/4) times their value: 
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With open boundary conditions this model shows SOC, 
while the Feder & Feder (FF) model which is iden- 
tical but for the increment, which is a constant that 
can be seen as the mean of aEi (pulse of the FF model 
= A = aEi), shows partial synchronization 11| (not 
SOC as claimed in This clearly means that the 

randomness of the initial conditions, which is dyamically 
eliminated in the FF model while it is mantained via the 
increment in the OFC model, changes the behavior of 
the system from partial synchronization to SOC. Further- 
more, different kinds of perturbations incompatible with 
a periodic behavior, change periodically ordered states 
for SOC. If, for instance, a random noise is added to the 
increment in the FF model, synchronization disappears 
and the system becomes SOC [0. On the other hand, if 
instead of open boundary conditions, periodic conditions 
are used, SOC disappears in favor of a periodic state, the 
period being the number of sites of the lattice It 
reappears if one inhomogeneity is introduced on the lat- 
tice |l^] . Extensive numerical simulations show that 
large avalanches are spatially quasi-stable. This means 
that during activation cycles of the oscillators, sites par- 
ticipating in large avalanches are synchronized. Thus the 
system finds a compromise between synchronization and 
SOC although this seems contradictory since SOC is in- 
compatible with complete synchronization Jl3| , pl| . 

In this letter I have shown that in the mean-field, re- 
laxation oscillator models such as the IF phase advance 
model synchronize under not as constrained conditions 
as suggested by the usual (erroneous) interpretation of 
MiroUo & Strogatz theorem. Moreover, disorder, which 
is a constant of the real world, does not spoil synchro- 
nization in this model, opening the prospect of a larger 
applicability of IF oscillators models and of variants. 
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FIG. 1. Plot of the probability Pjv, Eq.(6) for different dis- 
sipation levels. From top to bottom a — 0.5, 0.4, 0.3, 0.2, 0.1. 
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